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. Abstract 

There are major advantages in a newer version of Grover's quantum algorithm Q utilizing 
a general unitary transformation in the search of a single object in a large unsorted database. 
In this paper, we generalize this algorithm to multiobject search. We show the techniques to 
fT^ , achieve the reduction of the problem to one on an invariant subspace of dimension just equal to 

two. 
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1 Introduction 



This paper is a continuation from on quantum computing algorithms for multiobject 
search. 

L.K. Grover's first papers ^ on "quantum search for a needle in a haystack" have stim- 
ulated broad interest in the theoretical development of quantum computing algorithms. Let an 
unsorted database consist of N objects {wj | 1 < j < N}; each object Wj is stored in a quantum 
computer (QC) memory as an eigenstate \wj), j = 1,2, .. . ,N, with B = {\wj) \ 1 < j < N} 
forming an orthonormal basis of a Hilbert space Ti.. Let \w) be an element of B which is the 
(single) object to be searched. Grover's algorithm in H, ||] is to utilize a unitary operator 

U = -hlu, (1.1) 

where 



= / — 2|w)(w|, (J = the identity operator on "W) (1.2) 

1 ^ 

/,^7-2|.)(.|, (1.3) 

to perform the iterations C/™|s), which will lead to the target state \w) with probability close 
to 1 after approximately '^\/N number of iterations. The algorithm is of optimal order. 



In a more recent paper [Q, Grover showed that the state \s) in (1.3) can be replaced by any 
quantum state I7) with nonvanishing amplitude for each object Wj and, correspondingly, the 
Walsh-Hadamard operator previously used by him to construct the operator can be replaced 
by a sufficiently general nontrivial unitary operator. Grover's new "search engine" in Q is a 
unitary operator taking the form 

U = -I^V-^I^V: n^n (1.4) 

where V is an arbitrary unitary operator. The object w will be attained (with probability close 
to 1) by iterating ?7'"|7). 

This seems to give the algorithm/software designer large flexibility in conducting quantum 
computer search and code development. It increases the variety of quantum computational 
operations that can feasibly be performed by practical software. In particular, it opens the 
possibility of working with an initial state I7) (in place of \s)) that is other than a superposition 
of exactly N = 1"" {n = number of qubits) alternatives. This suggests a new paradigm in which 
the whole dataset (not just the key) is encoded in the quantum apparatus. This new point 
of view may also overcome some of the practical difficulties noted by Zalka in searching a 
physical database by Grover's method. 

In the next section, we study the generalization of (|l.4| ) to multiobject search. 
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2 Multiobject Search Algorithm Using a General Unitary 
Transformation 



Let {\wi) I 1 < i < N} be the basis of orthonormal eigenstates representing an unsorted 
database Wi, 1 < i < N, as noted in §1. We inherit much of the notation in 0]: let / be an 
oracle function such that 

^("^') = | 0, £+l<z<N, 
where Wi, i — 1,2, .. . ,£, represent the multiobjects under search. We wish to find at least one 
Wi, for i — 1,2, .. . ,£. Let I7) be any unit vector in Ti., and let L = span{|z«i) | 1 < i < ^}. 
Define 

j^=j-2|7)(7|: n^n, 

and 

/i|u;,>-(-l)^(-^)K), :j^l,2,... ,N, 

and II is then uniquely extended linearly to all Ti with the representation 

e 

/i = /-2^|u.,)(u;,|. 

Both I-y and II are unitary operators. Let V be any unitary operator on Ti.. Now, define 

U = -I-yV-^lLV. (2.1) 



Then J7 is a unitary operator; it degenerates into Grover's operator U in (1.4) when £ = 1 and 



further into the old Grover's operator U in (1.1) if y = / 



The unit vector I7) e is arbitrary except that we require V\^) ^ L. (Obviously, any I7) 



such that (w^Ij) ^ for alH = 1, 2, . . . ,N, wiU work, including I7) = \s) in (|l.3|).) If V\j) G L, 
then 

e I 

A measurement of the state VYi) will yield an eigenstate \wj), for some y. 1 < j < £, with 
probability \gj\'^. Thus the search task would have been completed. Thus, let us consider the 
nontrivial case V\'^) ^ L. This implies I7) ^ V^^{L) and, hence, 

Z = span({|7)}Ul^-i(L)) (2.2) 

is an {£ + l)-dimensional subspace of H. It effects a reduction to a lower dimensional invariant 
subspace for the operator U, according to the following. 

Lemma 2.1. Assume that (7I7) = 1 and V\j) ^ L. Then U{L) = L. 
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Proof. For any j : I < j < (, denote 



(1) We have, for j: I < j < i, 



U{V-'\w,)) ^ -L,V-' (^I-2j2H){w.,\^ \w,) 



(2) 



= -i^v-'Hw,)) 

= I^V-'\w,) 

= {I~2\^){^\)V-'\w,) 

^V-'\w,)-2i(j\V-'\w,)M 



U\^) = -I,V-' (^I-2j2\m){w..\j {V\j)) 

h)-2j2i{m\VmV'' 



4=1 



-a-2|7)(7l) 
e 

I7) +2^H^,,V-^\W^) -4.J2^^1,^'P^,^\l} 



Wi 



i=l 



(2.3) 



(2.4) 



By Lemma 2.1, the Hilbert space Ti admits an orthogonal direct sum decomposition 



such that L-^ is also an invariant subspace of U. In our subsequent iterations, the actions of 



U will be restricted to L, as the following Lemma 2.2 has shown. Therefore we can ignore the 
complementary summand space L-^. 



Lemma 2.2. Under the same assumptions as Lemma 2.1, we have U"^\'y) G L for m G 
{0,1,2,...}. 



Proof. It follows obviously from by (2.2) and Lemma 2.1 



□ 



Consider the action of U on L. Even though I7), V ^\wi), i = 1, . . . ,£, form a basis of L, 
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these vectors are not mutually orthogonal. We have 



U 



■ l7> ■ 
V-^\wi) 

V'^\W2) 



4 = 1 



M 



-2^7,1 

-2^7/ 

■ l7) ■ 

V-^\wi) 

V-^\W2) 








■ l7> ■ 
V-^\W2) 



(2.5) 



according to ( pT^ ) and (2.4). Therefore, with respect to the basis {I7), V~^\'Wi) | i = 1, . . . ,^}, 
the matrix representation of J7 on L is A^^, the transpose of M. These two + 1) x + 1) 
matrices and are nonunitary, however, because the basis {I7), y^^jwi) |, i = 1, 2, . . . , ^} 
is not orthogonal. This fact is relatively harmless here, as we can further effect a reduction of 
dimensionality by doing the following. Define a unit vector 

1/2 



\fi)=2^fi^^,V-'\w,)/c 



4l]lM7jl 



> 0. 



(2.6) 



Theorem 2.3. Let V = span{\j), \fJ.)}. Then V is a two-dimensional invariant subspace of U . 
We have 

1 — a 



U 



\l) 
1/^) 



= M 



l7) 



M = 



1 



(2.7) 



Consequently, with respect to the basis {I7), in V, the matrix representation of U is 



Proof. Using (2^), we have 

i e 
U\fi) = 2^/i7,,F->,) • - -2J2M' ■ -I7) 



3 = i 



Again, from the definition of in (2.6), we see that (2.4) gives 

U\J) = il-a'M+a\^i). 



Therefore (2.7) follows. 



□ 



Theorem 2^ gives a dramatic reduction of dimensionality to 2, i.e., the dimension of the 
invariant subspace V. Again, we note that the matrices M and AI^ in (2/7) are not unitary. 
Any vector \v) (z V can be represented as 



\v) = C1I7) +C2|Ai), 
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and so 



and thus 



C/|«) = C/(ci|7)+C2|m) 

= ci[(l - a^)\j) + a|M)] + C2[-a|7> + 



U\v) = A'f 



Cl 






—a 




Cl 


C2 




a 


1 




C2 



(2.8) 



where the first component of the vector on the riglit hand side of (2.8) corresponds to the 
coefficient of I7) while the second component corresponds to the coefficient of Therefore 

(2.9) 



"1 - 


—a 


m 


'1 


a 


1 








t/"|7) 

The above can be viewed geometrically (]5[) as follows 

[0,1] I ^> 




Figure 2.1: A geometric view of a single iteration ( p.§| ) 

M'^ [q] ~ o > very small, a « sin a, and therefore [^~°^] is a vector obtained 

from the unit vector [q] by rotating it counterclockwise with angle a. It takes approximately 

1 1/2 

e 

rotations to closely align the vector U™\'^) with e V^^L^. Thus V{U"''\j)) deviates little 
from the subspace L = spanjliUi) | i = 1, 2, . . . ,£}. A measurement of VU™\j) gives one of the 
eigenstatcs \wj), for some j: I < j < £, with probability nearly equal to 1, and the task of 
multiobject search is completed with this large probability. 
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